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Abstract 

We estimate the QCD effective charge ag in the low-energy region by exploiting 
I the conventional meson spectrum within a relativistic quantum-field model based on 

' analytic confinement. The ladder Bethe-Salpeter equation is solved for the masses 

^ . of two-quark bound states. We found a new, independent and specific infrared-finite 

^3^1 behavior of QCD coupling below energy scale 1 GeV. Particularly, an infrared-fixed 

^ ■ point is extracted at as{0) ~ 0.757 for confinement scale A = 345 MeV. As an appli- 

. cation, we estimate masses of some intermediate and heavy mesons and obtain results 

in reasonable agreement with recent experimental data. 
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1 Introduction 

The study of QCD behavior at large distances is an active field of research in particle 
OO '. physics because many interesting and novel behaviors are expected at low energies below 
1 GeV (see, e.g., [1, 2]). Understanding of a number of phenomena such as quark confine- 
^ ment, hadronization, the effective coupling and nonvanishing vacuum expectation values, 

Q ■ etc. requires a correct description of hadron dynamics in the infrared (IR) region. However, 
O . the well-established conventional perturbation theory cannot be used effectively in the IR 
region and it is required either to supply some additional phenomenological parameters (e.g., 
"effective masses", anomalous vacuum averages, etc.), or to use some nonperturbative meth- 
^ I ods (lattice simulations [3], power correction [4], string fragmentation [5], Dyson-Schwinger 
■ equations, etc.). There exists a phenomenological indication in favor of a smooth transition 
from short distance to long-distance physics [4]. 

One of the fundamental parameters of nature, the QCD effective coupling can provide 
a continuous interpolation between the asymptotical free state, where perturbation theory 
works well, and the hadronization regime, where nonperturbative techniques must be em- 
ployed. 

QCD predicts the functional form of the energy dependence of as on energy scale Q, 
but its actual value at a given Q must be obtained from experiment. This dependence is 
described theoretically by the renormalization group equations and measured at relatively 
high energies [6, 7]. A self-consistent and physically meaningful prediction of the QCD 
effective charge in the IR regime remains one of the actual problems in particle physics. 

The present paper is aimed to determine the QCD effective charge in the low-energy 
region by exploiting the hadron spectrum. In doing so we extend our previous investigations 
[8, 9, 10], where we provided new, independent, analytic and numerical estimates on the 
lowest glueball mass, conventional meson spectrum and the weak decay constants by using 
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a fixed ("frozen") value of as- The obtained results were in reasonable agreement with 
experimental evidence. 

Below we take into account the dependence of as on mass scale M and develop a phe- 
nomeno logical model to describe the IR behavior of as- We determine the meson masses by 
solving the ladder Bethe-Salpeter (BS) equations for two-quark bound states. The consid- 
eration is based on a relativistic quantum-field model with analytic confinement (AC) and 
has a minimal number of parameters, namely, the confinement scale A and the constituent 
quark masses m/, (/ = {wrf, s, c, 6}). First, we derive the meson mass formula and adjust 
the model parameters by fitting heavy meson masses (M > 2 GeV). Hereby, we determine 
corresponding values of as{M) from a smooth interpolation of the newest experimental data 
on the QCD coupling constant. Having adjusted model parameters, we estimate as(M) 
in the low-energy domain by exploiting meson masses below ~ 1 GeV. As an application, 
we estimate some intermediate and heavy meson masses (1 < M < 9.5 GeV). Finally, we 
extract a specific IR-finite behavior of the QCD coupling and conclude briefly recalling the 
comparison with often quoted results and recent experimental data. 

2 Effective Coupling of QCD 

The polarization of QCD vacuum causes two opposite effects: the color charge g is screened 
by the virtual quark-antiquark pairs and antiscreened by the polarization of virtual glu- 
ons. The competition of these effects results in a variation of the physical coupling under 
changes of distance ~ so QCD predicts a dependence as = g^/i^An) = as{Q)- This de- 
pendence is described theoretically by the renormalization group equations and determined 
experimentally at relatively high energies [6, 7]. 




Figure 1: Measurements of as as a function of the respective energy scale Q versa QCD 
predictions (curves) [11]. 

Nowadays, determinations of as remain at the forefront of experimental studies and tests 
of QCD. Recent developments on this subject were summarized in a number of articles 
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[2, 11, 12]. Summary of the recent experimental measurements of (Fig. 1) and particular 
values of as at intermediate energies (Tab. 1) are given by referring to [7, 11]. 

Note that there are two separate scale regions in which a running coupling may be 
considered. The spacelike region (g^ = — < with relativistic momentum transfer Q^) 
is related to scattering processes while timelike domain (g^ = > 0, where M is the 
hadron mass) is often used for annihilation and decay processes. The consistent description 
of QCD effective coupling ag in these domains remains the goal of many studies because 
only asymptotically the two definitions can be identified but at low momentum they can be 
very different (see, e.g. [13]). 

Particularly, the behavior of one-loop analytic running coupling [14] in timelike and 
spacelike domains is plotted in Fig. 2. 

Many quantities in hadron physics are affected by the IR behavior of the coupling in 
different amounts. Nevertheless, the long-distance behavior of as is not well defined, it 
needs to be more specified [15, 16, 17] and correct description of QCD effective coupling in 
the IR regime remains one of the actual problems in particle physics. Particularly, one of the 
most precise determinations of near the low-energy region is done by studying r-lepton 
decays reporting central values ranging from 0.318 to 0.344 [18, 19, 20]. 




Figure 2: The one- loop massless analytic running coupling in the spacelike and timelike 
domains (taken from [14]). 



An attempt to extrapolate the perturbative approach to the long-distance QCD has 
been made, it has been suggested that as freezes at a finite and moderate value [24], and 
this behavior could be the reason for the soft transition between short and long distance 
behaviors. 



Process 


Q [GeV] 


asiQ) 


Ref: 


T-decays 


1.78 


0.330 ± 0.014 


[11] 


QQ states 


4.1 


0.239 ± 0.012 


[21] 


T decays 


4.75 


0.217 ± 0.021 


[22] 


QQ states 


7.5 


0.1923 ± 0.0024 


[11] 


T decays 


9.46 


0.184 ± 0.015 


[11] 


e'^e" jets 


14.0 


0.170 ± 0.021 


[23] 



Table 1: Some measurements of as at intermediate energies. 
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Different nonperturbative approaclies fiave been proposed to deaf witli tlie IR properties 
of as- Particularly, methods, based on gauge-invariant SDE, concluded that an IR- finite 
coupling constant may be obtained from first principles [25] . New solutions for the gluon and 
ghost SDE have been obtained with better approximations which led to a new value for the 
IR coupling constant at the origin [26, 27]. Many works within the lattice simulations have 
been devoted in recent years to the study of the QCD running coupling constant either in the 
perturbative regime [28, 29] or in the deep IR domain [30]. Note that the results of various 
nonperturbative methods for the QCD invariant coupling may differ among themselves in the 
IR region due to the specifications of the used methods and approximations. Particularly, 
the results obtained by lattice simulations and SDE methods demonstrate a considerable 
variety of IR behaviors of as- 

An extraction of experimental data of al^^{Q'^) below 1 GeV compared with the meson 
spectrum within analytic perturbation theory has been performed [31] and a summary of 
data was presented (sec Fig. 2). The earliest attempts to obtain in the IR region were 
made in the framework of the quark- ant iquark potential models by using the Wilson loop 
method [32, 33, 34, 35, 36]. Convenient interpolation formulas between the large momentum 
perturbative expression and a finite IR- fixed point have been used in hadron spectrum studies 
with Q;°/7r ~ 0.19 — 0.25 [36]. Within a fully relativistic treatment it was shown that a p- 
meson mass much heavier than the tt mass could be obtained with a^/vr ^ 0.265 [37] while 
a similar result within a one-loop analytic coupling method predicted a^/vr ^ 0.44 [38]. A 
phenomenological hypothesis was adopted that the gluon acquires an effective dynamical 
mass rUg ^ 37QMeV (at Aqcd ~ SQQMeV) that resulted in a°/7r ~ 0.26 [39]. Various 
event shape in e'^e~ annihilation can be reproduced with an averaged value (cc^/tt) ~ 0.2 on 
interval < 1 GeV [4]. 

3 Model 

Color confinement in QCD is an attempt to explain the physics phenomenon that color 
charged particles are not observed. However, the reasons for quark confinement may be 
somewhat complicated. Particularly, within a quantum-field model, the quark confinement 
may be explained as the absence of quark poles and thresholds in Green's function. Following 
this idea, the conception of AC assumes that the QCD vacuum is realized by the self-dual 
vacuum gluon fields which arc stable versus local quantum fluctuations and related to the 
confinement and chiral symmetry breaking [40]. This vacuum gluon field serves as the true 
minimum of the QCD effective potential [41]. The vacuum of the quark- gluon system has 
the minimum at the nonzero self-dual homogenous background field with constant strength. 
Then, the quark and gluon propagators in the background gluon field represent entire analytic 
functions in Euclidean space [42]. In previous papers [43, 10]] we developed relativistic 
quantum-field models with AC. Similar ideas have been realized in infrared confinement by 
introducing an IR cutoff within a Nambu-Jona-Lasino model [44, 45]. 

The Bethe-Salpeter equation is an important tool for studying the relativistic two-particle 
bound states in a field theory framework [46]. Numerical calculations indicate that the ladder 
BS equation with a phenomenological model can give satisfactory results (for a review, see 
[47]). Particularly, a BS formalism adjusted for QCD was developed to extract values of ag 
below 1 GeV by comparison with known meson masses [31]. 

Our purpose is to investigate QCD effective (running) charge in the low-energy levels by 
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exploiting the spectrum of conventional mesons. For the spectra of two-quark bound states 
we consider a relativistic quantum-field model based on analytic (or infrared) confinement 
and solve the ladder BS equation. 

Following previous papers [43, 10] we consider a model Lagrangian: 



^--\ - gf^'^K^"".)" + E had- - + gr-c^^Y' gj) , (1) 



where A'^ is the gluon adjoint representation (a = {1, ...,4}); = d'^A^ — d'^A^; f'^^'" is 
the 5*^4(3) group structure constant ({A, -B, C} = {1, 8}); q'j is the quark spinor of flavor 
/ with color a = {1,2,3} and mass m/; g is the coupling strength, Fg = ija't''^] and t*" is 
the Gell-Mann matrices. 

Remember, that within the model the quark and gluon propagators S{p) and D(p) in 
(1) are entire analytic functions in the Euclidean space. 



3.1 Confinement and Green's Functions 

The effective charge is strongly governed by the detailed dynamics of the strong interaction 
and may depend on some of the most fundamental Greens functions of QCD, such as the 
gluon and quark propagators [48]. The Green's functions in QCD are tightly connected to 
confinement and are ingredients for hadron phenomenology. However, any widely accepted 
and rigorous analytic solutions to these propagators are still missing. One may encounter 
difficulties by defining the explicit quark and gluon propagator at the confinement scale. 
Nowadays, IR behaviors of the quark and gluon propagators are not well-established and 
need to be more specified [15]. 

The matrix elements of hadron processes at large distance are integrated characteristics 
of the vertices, quark and gluon propagators and the solution of the BS equation should 
not be too sensitive on the details of propagators. Taking into account the correct global 
symmetry properties and their breaking (also by introducing additional physical parameters) 
may be more important than the working out in detail of propagators (e.g., [49]). In previous 
papers we exploited simple forms of quark and gluon propagators [10, 43] which were entirely 
analytic functions in Euclidean space and behaved similarly to the explicit propagators 
dictated by AC [42]. 

Following [10] we introduce the quark propagator as follows: 



(2) 



where p = p^^^f^ and uj{z) = (1 + z /A) . The sign "±" corresponds to the self- and antiself- 
dual modes of the background gluon fields. In (2) chiral symmetry breaking is induced by 
AC. The interaction of the quark spin with the background gluon field results in a singular 
behavior S'±(p) ~ 1/m/ in the massless limit m/ — >■ 0. This expresses the zero- mode solution 
(the lowest Landau level) of the massless Dirac equation in the presence of an external gluon 
background field and generates a nontrivial quark condensate [10] indicating the broken 
chiral symmetry as m/ — )■ 0. 

Recent theoretical results predict an IR behavior of the gluon propagator. A gluon prop- 
agator identical to zero at the momentum origin was considered in [50, 51] while another 
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propagator was of order l/m^ [4], where nig is the dynamical gluon mass [52]. A renormal- 
ization group analysis [53] and numerical lattice studies simulating the gluon propagator are 
consistent with an IR- finite behavior [54] . We consider a gluon propagator 

^ 

It represents a modification of gluon propagator defined in [10] and exhibits an explicit 

IR-finite behavior D{0) ~ 1/A^. For simphcity D(p) in (3) is given in Feynman gauge. 
Note that the propagators in Eqs. (2) and (3) do not have any singularities in the finite 
plane in Euclidean space, thus indicating the absence of a single quark (gluon) in the 

asymptotic space of states. In fact, an IR parametrization is hidden in the confinement scale 

A. 



3.2 Two-quark Bound States 

We allow that the coupling remains weak (q;^ < 1) in the hadronization region. Then, the 
consideration may be restricted within the ladder approximation sufficient to estimate the 
meson spectrum with reasonable accuracy. The leading-order contribution to the two-quark 
(qq) bound states is determined by the partition function 

Z,^^ JjVqVqexpl^-{qS-'q) + ^{{qTAq){qTAq))^^ , 

{{•))d= fo^e-^(-^^-'-^)(.). (4) 



Our model has a minimal number of parameters, namely, the scale of confinement A and 
the constituent quark masses {rriud^^si^ci^h)- 

Below we briefly introduce the basic steps entering into our model on the example of the 
quark-antiquark bound states [10] defined by Zqq in (4). 

First, we allocate the one- gluon exchange between colored biquark currents 

£2 = y ^ jj dxidx2 (qf^{xi)i-f^t\^{xi)) D^^{xi,X2) {qf^{x2)nut^qf^{x2)) . (5) 
/1/2 

and isolate the color-singlet combinations. We perform a Fierz transformation 

(^7m)'^^'(^7.) = J2Cj-0j0j, J^{S, P, V, A, T} , 
J 

where Cj = {1,1,1/2,-1/2,0} and Oj = {/, ^75, i7^, 757^, ^[7^, 7^]/2}. For systems con- 
sisting of quarks with different masses it is important to go to the relative co-ordinates in 
the center-of-masses system and introduce the relative masses C,i = mi/ {mi + 7722). Then, 
introduce a system of orthonormalized basis functions {Uq{x)}, where Q = {rir, I, n} are the 
radial, orbital and magnetic quantum numbers. Diagonalize C2 on basis {Uq{x)} and use a 
Gaussian path-integral representation for the exponential 



6 



by introducing a colorless biquark current Jj^ and auxiliary meson fields Bj^f with H — 
{Q,JJij2}. Then 

= (^JJVqDqexp {-{qS-\) + g{BMJ^f)} 

By taking explicit path integration over quark variables we obtain 

^ Z = (exp {IVln [1 + ^(S^y^)^]})^ , 

where Vj\/{x, y) is a vertex function. 

Introduce a hadronization ansatz and this will identify Bji/{x) with meson fields carrying 
quantum numbers jV. We isolate all quadratic field configurations (~ B]^) in the "kinetic" 
term and rewrite the partition function for mesons [10]: 



N I ^ NAT' ) 



(6) 



where the interaction between mesons is described by the residual part Wres{BM\ ~ ^{B%)- 
The leading-order term of the polarization operator is 

as\MM'{z) = jjdxdy Um{x) asXjj'{z, x, y) Uj^>{y) , (7) 

and the Fourier transform of its kernel reads 

asXjj'{p, X, y) jdz e'^''Xjj'{z, x, y) 

^ l^V^y^(^)^ [OA, {k + 6p) OA, {k - (8) 

where Tr = TrcTr^ J2±'i Trc and Tr^ arc traces taken on color and spinor indices, correspond- 
ingly, while Y,± implies the sum over self-dual and antiself-dual modes. 
We diagonalize the polarization kernel on the orthonormal basis {Uj^}: 

jjdxdyU^;{x)\jj,{p,x,y)UM'{y) = 5^^' X^i-p') 
that is equivalent to the solution of the corresponding ladder BS equation. We rewrite 

Aaa(V) = ^Jd'k\Vj(k)fUj,(k,p), (9) 

Vj{k) = Jd''xUj{x)^/D{x)e-''^, 

^Ar{k,p) = ^Tr [OjSmi (k + iip) Oj,Sm^ (k - , 

where Vj{k) is a vertex and Il_^f{k,p) is the kernel of the polarization operator. 

In relativistic quantum-field theory a stable bound state of n massive particles shows up 
as a pole in the S-matrix with a center of mass energy. Accordingly, the physical mass of 
the meson may be derived from the equation 

l + aAM{M]^) = 0: -p^^Ml. (10) 
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Then, with a renormahzation 



dz 



the partition function takes the conventional form 

Z = jvBn exp |-i (Br [p" + M^-] Br) - Wres[BR]^ . 

3.3 Conventional Meson Spectrum and Running Coupling 



(12) 



We use the meson mass M as the appropriate characteristic parameter, so the couphng 
as{M) is defined in a timelike domain. On the other hand, most of known data on as{Q) 
are possible in the spacelike region. The continuation of the invariant charge from the 
spacelike to the timelike region (and vice versa) was elaborated by making use of the integral 
relationships between the QCD running coupling in Euclidean and Minkowskian domains 
(see, e.g. [55, 16]). 

Below we consider the most established sectors of hadron spectroscopy, the pseudoscalar 
P(0~"'") and vector V(l ) mesons. 

The dependence of meson masses on dig and other parameters is defined by Eq. (10). 
Note that the polarization kernel \m{~P^) is natively obtained real and symmetric that 
allows us to find a simple variational solution to this problem. Choosing a trial Gaussian 
function for the ground state [10] 



t/(x) = -^exp{-aAV} , A''Jd''x\U{x)f = 1, a>0. 
we obtain a variational form of Eq. (10) for meson masses as follows: 



-OisiMj) ■ \j{A,Mj, 1711,1712) 

8d,Cj f Mf - (mi + 1712? 



(13) 



(14) 



37r2(mi/A)(m2/A) 



• exp 



■ max 

0<c<2 



c (2 - c) 



2A2 

du dw 



' 2 p. J Mj 



Q 



A2 



, ^{l/u-l){l/w-l)Q 
66 + ^^^^^ 1 1 



■ exp 



(^1 - ^2 



2g 



2g. 



17111712 



4A2g j 



where Q ^ 1 + c{u + w), pj = {1, 1/2} and xj = {1, -1} for J = {P, V}. 

Further we exploit Eq. (14) in different ways, by solving either for otg at given masses, or 
for Mj at known values of coupling. In doing so, we adjust the model parameters by fitting 
available experimental data. 

Note that any physical observable must be independent of the particular scheme and mass 
by definition, but in (14) we obtain as depending on scaled masses {Mj/A, mi/A and m2/A}, 
where A is the scale of confinement. This kind of scale dependence is most pronounced in 
leading-order QCD and often used to test and specify uncertainties of theoretical calculations 
for physical observables. Conventionally, the central value of Q!s(/x) is determined or taken 
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for n equaling the typical energy of the underlying scattering reaction. There is no common 
agreement of how to fix the choice of scales. Particularly, in [10] we fixed the parameter A 
by fitting light meson weak decay constants. 

Below we solve Eq. (14) for different values of confinement scale. As a particular case, 
first we choose Ai = 345 MeV. 

1) We can extract intermediate values of as{My) in interval 2 — 10 GeV from a smooth 
interpolation of known data from Table 1. Particularly, 



a,(9460) = 0.1817, 
ds(3097) = 0.2619, 
0:5(2112) = 0.3074, 
0:5(2010) = 0.3138. 



(15) 



Hereafter, masses are given in units of MeV. 

Then, we adjust the constituent quark masses {'mud,'>Tis,'>Tic,mb} by solving a set of 
equations: 

1 + ^5(9460) ■ Ay(Ai, 9460, rrib, m^) = , 
1 + 0^(3097) ■ Ay(Ai, 3097, m^, m^) = , 
1 + 0:5(2112) • Ay(Ai, 2112, m„ rric) = , 
1 + 0:5(2010) • Ay(Ai, 2010, mud, rric) = 

with known masses of mesons T(9460), J/*(3097), D*(2112) and D*(2010). We fix a par- 
ticular set of model parameters as follows: 



(16) 



rus = 293.45 MeV , 



A = Ai = 345 MeV , 
mc = 1447.59 MeV , 



ruud = 192.56 MeV, 
rrib = 4692.51 MeV . 



(17) 



2) Having fixed the model parameters, we solve an inverse problem, to find 0:5 values in 
the region below 1 GeV as follows: 



65(138) = -\p\Ai, 138, mud, mud) = 0.7131 , 
0:5(495) = -Ap^(Ai,495,m„d,TOs) = 0.6086, 
65(770) = -Xy\A^, 770, mud, mud) = 0.4390, 
[ 65(892) = -A^^(Ai,892,mw,m5) = 0.4214. 



(18) 



In Fig. 3 we plot our low-energy estimates (18) in comparison with the three-loop ana- 
lytic coupling, its perturbative counterpart (both normalized at the Z-boson mass), and the 
massive one-loop analytic coupling [31]. 

3) As an application, with particular choice of parameters (17) we calculate masses of 
other mesons: D(1870), L>5(1970), 77c(2980), 5(5279), S*(5325), ^5(5370), Se(6286) and 
77^(9389). Hereby, the corresponding 0:5 (M) are extracted from Fig. 1. 

Our estimates of meson masses along experimental data [2] are shown in Table 2. The 
relative error of our estimate does not exceed 3.5% percent in a wide range of mass. 

4) To check the sensibility of the obtained results on the confinement scale value we 
recalculated steps 1-3 for A = 330 MeV and A = 360 MeV. We revealed that the estimated 
meson masses shown in Table 2 do not change considerably (less than 0.5% percent). The 
variation of 65 under changes of A is shown in Fig. 3. 

5) Wc perform global evaluation of 65 (M) at the mass scale of conventional mesons 
(shown in Table 2) by using the formula 



65(Mj) = -l/Aj(Mj,A,mi,m2) 
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Figure 3: Our estimates of as{M) in the low-energy region at different values of confinement 
scale (red dots for A = 330 MeV; blue diamonds for A = 345 MeV and black squares for 
A = 360 MeV) compared with the three- loop analytic coupling as{Q) (solid curve), its 
perturbative counterpart (dot-dashed curve), and the massive one- loop analytic coupling 
(dashed curve) (see Ref. [31]). 

and we plot the resulting curves at different A in Fig. 5 in comparison with recent low- and 
high-energy data of as{Q) [31]. 

3.4 IR-finite Behavior of Effective Coupling 

The possibility that the QCD coupling constant features an IR-finite behavior has been 
extensively studied in recent years (e.g., [56, 57]). There are theoretical arguments in favor 
of a nontrivial IR- fixed point, particularly, the analytical coupling freezes at the value of 
4:71 /Po within one- loop approximation [58]. The phenomenological evidence for as finite in 
the IR region is much more numerous. 

We note that the agreement of our estimates of as(M) with other predictions (e.g., [7, 13]) 
turns out to be reasonable from 2 GeV down to the 1 GeV scale. Below this scale, different 
behaviors of as{M) may be expected as M approaches zero. 

Below we consider the IR-fixed point d° = ds(0) by evaluating Eq. (14) for Mp = and 
mi = 7712 = rn: 



~8A^' 



max 

0<c<2 



c(2 



du dw 



i ^(l/u-l)il/w-l)il + ciu + w)y 



{i + c{u + w)y 



(19) 



The dependence of d° on = m/A is plotted in Fig. 4. 

Note that a value of of order 2 or larger would be definitely out of line with many other 
phenomena, such as nonrelativistic potentials for charmonium [59] and analytic perturbation 
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J^^' = 0-+ 


Mp 


J^^' = 0-+ 


A'/p 


J^^' = 1- 


Mv 


J^^' = 1 


Mv 


7r(138) 


138 


r/c(2980) 


3039 


p(770) 


770 


L);(2112) 


2112 


K(495) 


495 


fi(5279) 


5339 


w(782) 


785 


J/^(3097) 


3097 


r/(547) 


547 


5,(5370) 


5439 


K*(892) 


892 


5*(5325) 


5357 


D(1870) 


1941 


5c(6286) 


6489 


$(1019) 


1022 


T(9460) 


9460 


Dsimo) 


2039 


r/6(9389) 


9442 


L)*(2010) 


2010 







Table 2: Masses M of conventional mesons (in units of MeV) corresponding to effective 
coupling ds(M) determined by Eq. (14) at A = 345 MeV. 
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Figure 4: Dependence of IR fixed point a° on the scaled quark mass /i^ = (m/A)^ at fixed 
confinement scale A = 345 MeV. 



theory [58]. Obviously, this constraint implies an upper limit to the value of constituent quark 
mass: /i^ < 0.8 or m < 0.9 A. 

Since we are searching the IR- fixed point, it is reasonable to choose the lightest quark 
mass. Particularly, for m = iriud = 192.56 MeV and A = 345 MeV we obtain 

= 0.757 , or d^/vr = 0.241 . (20) 

To compare our result with known data on as{Q) we exploit the integral relationships 
between the QCD running coupling in Euclidean and Minkowskian domains. Particularly, 
there exists a relation [16] 

^si<l') = <l'j7^f-^,^s{s) (21) 
valid for the case of massless pion. By substituting s = t(f into (21) one rewrites 

oo 

«^(^')= / TTXT^"^^^^')- ^22) 
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Then, for — )■ we obtain 



a,(0) = d,(0) 



dt 



«s(0)-l. 



(23) 



(24) 



Therefore, we may conclude that our result (20) is in reasonable agreement with often- 
quoted estimates 

f a°/7r ~ 0.19 - 0.25 [36], 
ay 71 ~ 0.265 [37] , 

a°/7r ~ 0.26 [39] , 

(«°A)iGey^0.2 [4] 

and phenomenological evidences [38, 31]. The obtained IR-fixed value of the coupling con- 
stant is moderate, it depends on the mass of constituent quark (m, d), so one can insert this 
value into perturbative expressions to be compatible with the experimental data. 




Figure 5: Summary of estimates of as{M) in interval from to 10 GeV at different values of 
confinement scale. In the left panel A = 330 MeV (red dots), A = 345 MeV (blue diamonds) 
and A = 360 MeV (black squares) compared with as{Q) (in the right panel) defined in low- 
energy (open diamonds) and high-energy (open circles) experiments. Also shown are the 
three-loop analytic coupling (solid curve), its perturbative counterpart (dot-dashed curve) 
both normalized at the Z-boson mass, and the massive one-loop analytic coupling (dashed 
curve) (for details see Ref. [31]). 

By interpolating smoothly obtained results in (20), (18) and (15) into the intermediate- 
energy region we define on a wide interval 0.14 — 9.5 GeV. Some particular cases of the 
dependence as on mass scale M at different model parameters are plotted in Fig. 5. 

It is important to stress that we do not aim to obtain the behavior of the coupling 
constant at all scales. At moderate = —p^ we obtain as in coincidence with the QCD 
predictions. However, at large mass scale (above 10 GeV) as decreases much faster than 
expected by QCD prediction. The reason is the use of confined propagators in the form 
of entire functions, Eqs. (2) and (3). Then, the convolution of entire functions leads to a 
rapid decreasing (or a rapid growth in Minkowski space) of physical matrix elements once 
the hadron masses and energies of the reaction have been fixed. Consequently, the numerical 
results become sensitive to changes of model parameters at large masses and energies. 
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4 Conclusion 



To conclude, we provide an estimate of QCD effective charge in the low-energy region (below 
1 GeV) by exploiting the conventional meson spectrum within a relativistic quantum-field 
model based on analytic (or infrared) confinement. The new results obtained in the previous 
section arc summarized in Figs. 3-5 and Table 2. 

We demonstrate that global properties of the low-energy phenomena such as QCD run- 
ning coupling and conventional meson spectrum may be explained reasonably in the frame- 
work of a simple relativistic quantum-field model of quark-gluon interaction based on analytic 
(or, infrared) confinement. Our guess about the symmetry structure of the qTiark-gluon inter- 
action in the confinement region has been tested and the use of simple forms of propagators 
has resulted in quantitatively reasonable estimates. 

Despite its pure model origin, the approximations used, and questions about the very 
definition of the coupling in the IR region, our approach demonstrates a new, independent 
and specific IR-finite behavior of QCD coupling and we extract a particular IR-fixed point at 
ds(0) ~ 0.757 for confinement scale A = 345 MeV. As an application, we performed estimates 
on intermediate and heavy meson masses and the result was in reasonable agreement with 
experimental data. Our estimates may be improved further by using iterative schemes, but 
the aim is to obtain a qualitative understanding of QCD effective coupling in the IR region. 

The suggested model in its simple form is far from real QCD but we conclude that the 
analytic confinement conception combined with BS method may provide us with a rather 
satisfactory correlated understanding of low and intermediate-energy phenomena from few 
hundreds MeV to few GeV. 

Note that further improvements of measurements of will be difficult while it is un- 
likely that QCD perturbation theory will considerably improve existing predictions. There- 
fore, further developments of theoretical predictions within nonperturbative methods and 
reapplication of improved models may have successes in this field. 

The author thanks M.A. Ivanov, E. Klempt and A.V. Nesterenko for useful discussions 
and valuable remarks. 
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